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[57] ABSTRACT 

A method for adaptation to the time-variable channel 
impulse response during a burst in receivers for TDMA 
mobile radio systems, each burst including in the center a 
test sequence which is known to the receiver and by means 
of which an initial estimate of the channel impulse response 
is possible, wherein 

the received signal of a burst is suitably sampled, 

the received burst is divided in 2N+1 sectors A_^, . . . , A^ 
it being assumed that the channel impulse response can be 
regarded as being constant within a sector, 

a vector x of dimension (2N+1)L is denned, L being the 
number of channel impulse response coefficients and the 
(iLfl)-th component of x, -N^i^+N, 0^1^L-1 is initial- 
ized before the start of the detection using the 1-th coefficient 
of the estimated value of the channel impulse response 
obtained from the test sequence, Matrices E, E ( , E w S and 
vectors i„ d, r ( are defined, 

each burst is processed in 2N steps, where 

in the k-th step, the i-th sector is detected using the 
components iL, . . . ,(i+l)L-l of the vector x, where the 
sector A ( - is selected in accordance with a suitable 
method, for example that described in [5], 

the matrix 



is calculated, 

and the recursions 



Er=E+PSET 

are carried out. 

1 Claim, 3 Drawing Sheets 
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KALMAN FILTER FOR CHANNEL 
IMPULSE-RESPONSE ADAPTION IN 
RECEIVERS FOR TDMA MOBILE RADIO 
SYSTEMS 



A measure of the rate of change of the channel is the 
so-called Doppler frequency 



15 



20 



25 



30 



1. BACKGROUND OF THE INVENTION 

The invention relates to the problem of adaptive data 
detection in a mobile radio receiver. "Adaptive" means that 10 
a device is provided in the receiver, which device can 
understand the time changes of the properties of the mobile 
radio channel. Such devices normally operate using the 
LMS (Least Mean Squares) (see for example, citation [1]), 
or the RLS (Recursive Least Squares) (see for example, 
citation [2]) algorithm. These algorithms adapt to the prop- 
erties of the transmission channel continually with time. 

EP-A-520 969 specifies an estimation method for data 
detection, in the case of which received samples in a time 
slot are initially stored and processed starting with the first 
received value, to be precise to beyond a fading incursion; 
further processing is men carried out, starting with the last 
stored sample, in order to produce estimated values of the 
transmitted data. 

A Kalman filter for adaptation of the channel impulse 
response in a mobile radio receiver is described in Fifth 
European Signal Processing Conference, Sep. 18-21 1990, 
Barcelona, ES; Elsevier, Amsterdam, NL, 1990; pages 
225-228, Morgul & Dzung: "New optimum recursive 
parameter estimation/detection using unreliable erasure 
declaring detectors". In this case, the Kalman algorithm is 
applied only to the channel impulse response which is 
present at a specific time. 

The present patent application shows that the algorithms 35 
which are continuously adapted with time fail under certain 
conditions, which are typical for mobile radios. The reason 
for this is that such adaptation algorithms can be regarded as 
time-variant first order recursive filters. Higher order algo- 
rithms exhibit a behavior which is better in principle. 

FIG. 1 shows the data structure which are transmitted in 
TDMA bursts in a GSM mobile radio system. 

FIG. 2 shows a block diagram of the transmission path. 

FIG. 3 illustrates the splitting of a TDMA burst into 
sectors. 

FIG. 4 shows a fading incursion in a TDMA burst. 
FIG. 5 illustrates a process model relating to the Kalman 
theory. 

50 

2.1 Model of the transmission system 

In digital TDMA mobile radio systems, the data are 
transmitted in bursts. In the GSM system, 1 16 data bits, for 55 
example, are combined to form a burst. The burst structure 
which is chosen for this purpose is illustrated in FIG. 1. In 
addition to the 116 data bits, which are grouped into two 
sections of 58 bits each, the burst also includes at each of the 
ends 3 so-called "tail bits" and in the center a so-called test e 0 
sequence, comprising 26 bits, so that the burst comprises a 
total of 148 bits. 

The test sequence is used to "survey" the transmission 
channel. This transmission channel can be regarded as a 
linear time-variant filter [4] and surveying of the transmis- 65 
sion channel thus has the same meaning as an estimate of the 
channel impulse response. 



40 



45 



where f 0 is the carrier frequency of the transmitted signal, 
v the speed of the mobile subscriber and c the speed of light. 
The length of a burst is selected such that, despite the 
possible changes, the channel impulse response during a 
burst can be regarded as being virtually constant. In the case 
of GSM, this assumption is valid for speeds v of up to about 
250 km/h (the carrier frequency f 0 is around 900 MHz, and 
the burst duration is around 546 us). 

The method according to the invention is predicated on a 
burst transmission similar to that in the case of GSM. In 
particular, the presence of a test sequence in the center of the 
burst is significant. 

Those parts of the transmitter and receiver which are 
important for the definition of the signals are illustrated in 
FIG. 2. The modulator maps the data sequence a to be 
transmitted into the signal s(t). It is assumed that s(t) is a 
complex low-pass signal having the quadrature components 
s/t) and s G (t), that is to say 

j(0=*X'}+->c(0 

In this case, j=^-T is the imaginary unit It is furthermore 
assumed that is the bandwidth of the signal. s(t) can thus 
be represented unambiguously, in accordance with the sam- 
pling theorem, by the samples s(kT r ) where T=l/W/s(t) is 
mapped by the quadrature modulator into the RF signal 

where f 0 is the desired carrier frequency. s fl (t) is further 
amplified if necessary and transmitted via the antenna. 

r fl (t) is assumed to be the RF received signal to which, as 
s B (t), in accordance with 



ra(^{K'V^}=r//)cos2^bf-r c (0sin Inftf 

a complex low-pass signal 



r<0=r/'>+-;> c (0 

can be assigned. The quadrature demodulator at the receiver 
end is used to obtain this signal. Like s(t), r(t) is a low-pass 
signal of bandwidth W,; r(t) can thus likewise be represented 
by the samples r(kTJ. 

In mobile radio, it is assumed that the relationship 
between the received signal samples and the transmitted 
signal samples is: 



The sum describes a time-discrete convolution; h(lJcTs) is 
the 1-th coefScient of the (time-variant) channel impulse 
response at the time t=kT J ni(kT J ) symbolizes an additional 
interference component; in mobile radio, this is normally 
composed of thermal noise and signals from other transmit- 
ters. It is assumed in the following text that this interference 
component has no mean value, has a gaussian distribution 
and has the correlation properties 



06/25/2004; EAST Version: 1.4.1 



5,615,208 



<n(kT,)n*{lT t )>=a*SkL 

<x> is in general used as the notation for the expected value 
of the random variables x. 

The time constant is always omitted from the notation 5 
in the following analyses; it is furthermore assumed that the 
channel impulse response coefficients differ from 0 only for 
1=0, . . .L-l. The above relationship between the received 
signal samples and the transmitted signal samples then 



becomes 



10 



1=0 



(1) 



In order to simplify the notation, the channel impulse 
response coefficients at the time kTj are combined to form 15 
the vector 



30 



35 



h{k)=(h{Q t k), . . . ML-U)) T 

With respect to the notation: multidimensional variables 
such as vectors and matrices are marked by an underscore. 
The first summand in equation (1) can be called the useful 
component, and the second the interference component of 
the received signal. 

2.2 General principle of detection 

In order to describe the method according to the inven- 
tion, it is important that the channel impulse response in the 
receiver must be known in order to determine the transmit- 
ted data. This follows directly from equation (1). If it is 
assumed that the noise signal is not present, that is to say 
n(k)=0, then the received signal is dependent only on the 
transmitted signal and the channel impulse response. Thus, 
if the received signal and the channel impulse response are 
known, then equation (1) can be solved on the basis of the 
transmitted signal and the transmitted data can subsequently 
be determined by inversion of the modulation. The converse 
applies: if the data are known, then (1) can be solved on the 
basis of the channel impulse response. If n(k)*0, then it is 40 
not possible to find any direct solution of equation (1) on the 
basis of s(k-l) or h(l,k); nevertheless, it is also true here that 
the channel impulse response must be known in order to 
determine the transmitted signal, and vice versa. The deter- 
mined values for the data and the channel impulse response 45 
are estimated values, because the noise n(k) is not insignifi- 
cant. 

The general principle of detection can now be described 
as follows. The received burst is divided into 2N+1 sectors 
A,-, -N^i^+N, which do not necessarily need to be of the 
same length. However, the length of the sectors should be 
selected such that the channel impulse response can be 
assumed to be constant within a sector. The sector A 0 should 
be identical to the test sequence. Such splitting into sectors 55 
is illustrated in FIG. 3. 

Since the test sequence is known, the channel impulse 
response can be determined for the sector Aq. Since the 
channel impulse responses in adjacent sectors differ from 
one another only slightly, the data in the sectors A^ and A+i ^ 
can be determined with the aid of the estimated value of the 
channel impulse response in the sector Aq. The channel 
impulse responses of these sectors can now be estimated 
using these data. The channel impulse responses are then in 
turn used for data detection in the sectors A_ 2 and A^ 2 , etc. $5 

The estimation of the channel impulse response in the 
sector A ; , i>0 will be based not only on the data detected in 



this sector but also, for example, on the channel impulse 
response in the sector A^. The reason for this is that the 
channel impulse responses vary only slightly from sector to 
sector; the channel impulse responses of the individual 
sectors are highly correlated with one another. 

The idea on which the invention is based consists essen- 
tially of the knowledge that the correlations of the channel 
impulse response coefficients should be used as far as 
possible, that is to say all the channel impulse responses of 
the already processed sectors should be considered in the 
detennination of the channel impulse response for the sector 
A;. However, this makes sense only if the estimated values 
of the channel impulse responses of the already processed 
sectors are also reliable: For this reason, the processing of 
the burst should be carried out, starting with the most 
reliable sector, in the direction of reducing reliability (see 
[5]). 

2,3 Previously known solutions 

20 

The mostly widely known methods for channel impulse 
response adaptation are the LMS algorithm and the RLS 
algorithm. In both cases, the sectors outside the test 
sequence each consist of only one signal sample. h(i) is 
25 assumed to be the estimated value of the channel impulse 
response for the i-th sector. For i>0, both algorithms calcu- 
late 



50 



h(i>h(i-i)4-/>(0rf(i). (2) 
Where P(i) is a vector of the same dimension as h(i) and 



(3) 



is the error between the actual received signal r(i) and the 
received signal model 



L-l « 
K0= £ h(l,i- 



lMi-J). 



(4) 



Equation (2) describes the fact that h(i) is dependent only on 
n(i-l), that is to say the estimated channel impulse response 
in the preceding sector. The detected data are required to 
calculate the error d(i). 

These analyses apply both for the LMS algorithm and for 
the RLS algorithm. The difference between the two algo- 
rithms is the choice of P(i). Since every new estimated value 
of the channel impulse response depends, in accordance 
with equation (2), only on the immediately preceding value, 
it is possible to speak of first order algorithms (by analogy 
with the theory of recursive time-discrete filters). 

The time direction in the adaptation equation (2) can be 
reversed. This then gives 



h(0=h(£+i>fP(iV(0. 



(5) 



This is important for processing the sectors A, where i<0. 
First order algorithms fail under two conditions. 
1 . If the estimated value of the channel impulse response 
in sector A^ is completely false, then the data detec- 
tion in the. sector A f is also erroneous since this is 
actually being carried out using the estimated channel 
impulse response of sector A*.!. Since the estimated 
value of the channel impulse response in the sector A f 
depends only on h(i-l) and the data of sector A f , this 
is also completely false. An incorrect estimate of the 
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channel impulse response in the sector A,-^ thus leads, 
with high probability, to an incorrect estimate of the 
channel impulse response in the sector A,-. 
2. If the components of the channel impulse response 
vector in the sector A, are so small that the useful 5 
component of the received signal is small in compari- 
son with the interference component, then there is a 
high probability that the data detection in this sector 
will be erroneous. It now follows from equations (2), 
(3) and (4) that the change in the channel impulse 10 
response from the sector A,-.! to the sector A, cannot be 
understood completely; the estimated change is, rather, 
completely random. It follows from this and from the 
small magnitude of the components of h(i) that there is 
a high probability that the estimated value h(i) will be 15 
completely false. 
. From these two items, it follows that the first order 
algorithms do not permit adaptation resulting from a fading 
incursion. A fading incursion is understood to be a situation 
in which the channel impulse response coefficients assume 20 
such small values that the useful component of the received 
signal disappears in the noise. Let us assume that the sector 
Aj_ x has been affected by a fading incursion. It now follows 
from the second item that there is a high probability that the 
estimated value of the channel impulse response in the 25 
sector A,, will be highly erroneous. It now follows from the 
first item that this also applies to all further sectors A i+ll A^ 2 
etc. 

If a fading incursion occurs in a data section, then the data 
behind (seen from the test sequence) this fading incursion 30 
can no longer be detected. The probability of such a fading 
incursion rises with the vehicle speed. This is the reason why 
first order algorithms fail at high vehicle speeds. In the GSM 
system, this statement applies to vehicle speeds above 
approximately 250 km/h. 35 



which cannot adapt as a result of this fading incursion, the 
data bits at the end of the second data section would no 
longer be detectable. However, since the channel impulse 
response coefficients on both sides of the fading incursion 
are correlated with one another, it is quite possible to use the 
estimated channel impulse responses on the left-hand side to 
deduce those on the right-hand side of the fading incursion. 
However, it should be ensured for this purpose that the 
channel impulse response profile on the left-hand side of the 
fading incursion is completely known before any attempt is 
made to determine that on the right-hand side. A device must 
therefore be provided which initially provides all the reliable 
data regions, that is to say data regions which are not 
affected by a fading incursion, for detection before an 
affected region is processed (see equation [5]). 

Before the description of the actual method, a brief 
summary of the theory of the Kalman filter should be 
permitted. The process model on which the Kalman theory 
is based is a vectorial sequential state mode] (see FIG. 5). 
x(k) is the state vector, u(k) the control vector and y(k) the 
output vector for the recursion step k in each case. The 
process model states that the state vector in the recursion 
step k+1 depends linearly on the state vector and the control 
vector in the recursion step k, i.e. 



2.4 Channel impulse response adaptation with the 
aid of a Kalman filter 



40 



The possibility for adaptation of the channel impulse 
response resulting from a fading incursion is a precondition 
for reliable mobile communication at very high vehicle 
speeds. Examples of the use of this are fast-moving trains, 
for which speeds up to 500 km/h are being discussed. The 45 
GSM system was initially designed for vehicle speeds up to 
around 250 km/h. Use of the GSM standard for mobile 
railway radio is possible only if the methods for channel 
impulse response adaptation and, in particular, its capability 
to adapt the channel impulse response caused by a fading 50 
incursion can be significantly improved throughout. 

A significant improvement in the channel impulse 
response adaptation is possible only by means of higher 
order algorithms. In the case of such algorithms, the esti- 
mated value h(k) is dependent not only on h(k-l) (or h(k+l) 55 
in the case of backward adaptation), but on as many already 
determined values of h(k+n), r&Q as possible. This enables 
improved utilization of the correlations of the channel 
impulse responses at different times. 

Example: Let us assume the situation illustrated in FIG. 4. 60 



E h (t) 



L-\ 1 

/=0 



is the energy of the channel impulse response at time t. A 65 
fading incursion can be seen approximately in the center of 
the second data section. If a first order algorithm is used, 



(6) 



The linear mapping is described by the matrices A(k) and 
B(k) which, in general, depend on the recursion step. The 
output vector is linearly dependent on the state vector 



y(k)=C(k)x(k). 

The mapping matrix C(k) likewise in general depends on 
the recursion step. The initial state x(0) must be known in 
order to describe the recursion for k^0 completely. 

The object of the Kalman filter is to determine an esti- 
mated value for the state vector x(n+l) in the n-th recursion 
step, only disturbed versions of the output vectors y(0), . . , 
,y(n) being known to the filter. These disturbed versions are 
designated by z(k) and it is assumed to be a precondition that 
they are additively composed of the output vectors y(k) and 
disturbance vectors v(k), that is to say 



z(Jfc>=y(JtHv(Jt)=0(W*Hv(fc). 



a) 



The model is completely described by the equations (6) and 
(7). 

The theory of Kalman filters is based on the assumptions 
that the control vectors and disturbance vectors are random 
processes with no mean values, the vectors being uncorre- 
lated in different recursion steps, that is to say 



<u{k)>=<v(k)>=0 

<«(*)« *(D>«U(*)B„ 

<v(*)v*<o>=W8 u . 

The superscript * in this case means conjugation and trans- 
position of the designated vector (or matrix). The matrices 
U(k) and V(k) describe the correlations of the components 
in the vectors u(k) and v(k) £(k) is assumed to be the 
estimated value for x(k) for the (k-l)-th recursion step. 
Furthermore, it is assumed that 
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20 



is the autocorrelation matrix of the estimated error. The 
Kalman filter can now be described by the following recur- 
sive algorithm. 



p(kw{k)E(k)c*(k)mk}+c(k)E(k)c*m' 1 (8) 

i(k+\MMkyP(k)CXk))W)+P(k)z{k) (9) 

E{k+lMA(khP(k)CmmiMk)-P{k)C(k)}* +i><*)tW*(*)+ 10 
B{kmk)B*{k). (10) 

The first of these equations is not a recursion equation and 
only the auxiliary variable P(k) is calculated using it The 
actual recursion relates to the estimated value 5t(k) for the 
state vector x(k) and the error correlation matrix E(k). 
Suitable start values x(0) and E(0) must therefore be found 
for these two variables. 

The method according to the invention now results from 
the application of these three equations to a suitable model 
for channel impulse response adaptation. As is described in 
Section 2.2, the burst is split into 2N+1 sectors A r , i=-N, . 
. . ,+N, the assumed precondition here being that all the 
sectors are of the same length, that is to say each comprise ^ 
M signal samples. It is assumed that the sector A 0 is located 
completely in the test sequence, the channel impulse 
response for this sector is therefore known with the estima- 
tion of the channel impulse response from the test sequence. 
The transmitted signal samples of the i-th sector are s(iM), 
. . . ,s((i+l)M-l), the associated received signal samples are 
r(iM), . . . ,r((i+l)M-l), and the noise samples are n(iM), . 
. . ,n((i-f-l)M-l). These samples can be combined to form the 
vectors s ( , r, and n r -. Furthermore, it should be assumed that 
the channel impulse response is constant within a sector. 
This assumption becomes the more accurate the smaller are 
the sectors, that is to say the smaller M, the number of signal 
samples per sector, is. The channel impulse response vector 
of the i-th sector is now designated h,-, i=-N, . . . ,+N. 

The general sequence of processing a burst is character- 
ized by the following points. 

A burst is processed in 2N steps, the sector A ( which is 
provided for this purpose by the reliability-controlled 45 
data detection being detected in the k-th step [5]. 
An estimated channel impulse response h,(k) is assigned 
to each sector A t in the k-th step; this estimated value 
of the channel impulse response is used for the detec- 
tion of the sector. With regard to the notation: the index 
i designates the sector to be processed, the variable k in 
brackets designates the processing step. h,(k) is thus the 
estimated channel impulse response for the i-th sector, 
as results in the k-th processing step. The 1-th compo- 
nent of a vector fi is designated x(l). Theoretically, this 
can result in it being possible to confuse the processing 
step and the component. However, the meaning in each 
case becomes clear from the context. 
After the detection of the sector A f in the k-th step, the 
transmitted signal samples of this sector can be recon- 
structed; to put this more precisely, the detection must 
be carried out such that the transmitted signal samples 
s(iM-L+l), . . . ,s((i+l)M-l) can be reconstructed. The 
matrix of dimension MxL can now be produced using 
these samples: 



55 



60 



W(i + lW-l) 



s{iM-L+l) 



j((i4-l)M-L) 



(U) 



As will be shown, this matrix can now be used to calculate 
a recursion step of the Kalman filter which updates the 
channel impulse responses for all the sectors. 

An initial value for the channel impulse response is 
required not only for the Kalman recursion but also for 
the detection of the first sector. This initial value, which 
below will be designated h, can be obtained with the aid 
of the test sequence provided in each burst. This 
estimated value is initially assigned to all the sectors; 
that is to say h r -(0)=h applies for all i. 



(12) 



30 



35 



40 



is now defined as the state vector, that is to say the 
concatenation of the channel impulse responses to be esti- 
mated for all the sectors. This state vector does not, of 
course, change with k. Therefore, it immediately follows 
from equation (6) 



A(kH 
B(k)=0 
1100=0, 

where I is the unit matrix. All the variables in the first model 
equation (6) have therefore been found. The variables in the 
second model equation (7) are found by inserting equation 
into equation (1) (11) in; 



is obtained for the i-th sector. 
If one now defines: 

m = r, 

and 



(13) 



50 




then equation (7) is obtained from equation (13). Tne model 
is therefore completely described. It should be noted at this 
point that the correlation matrix of the model noise V(k) is 
equal to o^I since the model noise, as seen above, is equal 
to the channel noise but this has been assumed to be 
uncorrected with the variance a 2 . 

In addition to the model, the initial values x(0), £(0) and 
E(0) must be defined. x(0) results in a trivial manner from 
(12). £(0) is obtained, as described above, by assigning to all 
the sectors the estimated value of the channel impulse 
response h from the test sequence. 
Thus 



65 



*(0Mb T , 
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The initial value for the error correlation matrix can now be 
calculated from x(0) and x(0). 



£(0)=<[X(0>-Jt(0)]U(0)-x(0)]*>=<x(0)X*(0)>^0)x* (0)^- 
<*(0)X*(0>-cx(0)x* (0)>. 

This matrix has the dimension (2N+l)Lx(2N+l)L, its ele- 
ments being annotated in the following text as *Li+l,Li'+ 
r(0), -Ngyg+N, OlJ'SL-l. Thus, 



10 



^ +( .ii- + r(0)^(0&^/)>^^<0V"(O^,<0fi*(n>^(0Ar*(/>. 

If the test sequences have been suitably selected, then 

where n,,(l) is a random variable which is independent of 
h 0 (l) f has a gaussian distribution and has no mean, where 

<n A (JK*(/>=a A 2 5„. 

In practical applications, this assumption is valid at least to 
a good approximation. It follows that: 
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0. 



t» Wl (U'.i,0H» AA (U'.0, 



where 



40 



45 



10 

i the correlation of the 1-th component of the channel 
impulse response in the sector A, and the l'-th compo- 
nent of the channel impulse response in the sector A f , 
and o 2 /, is the variance of the error in the estimation of 
the channel impulse response coefficients from the test 
signal. E, is assumed to be a matrix of dimension 
(2N+l)LxL whose columns, after the detection of the 

i-th sector, are the columns iL (i+l)L-l of the 

matrix E. Furthermore, E„ is assumed to a matrix of 
dimension LxL whose lines, after the detection of the 
i-th sector, are the lines iL, . . .,(i+l)L-l of the matrix 

E* 

i matrix S of dimension MxL is defined, whose (mj)-th 
digit, after the detection of the i-th sector, is occupied 
by the reconstructed transmitted signal sample s(iM+ 
m-1). 

i vector % { of dimension L is defined, whose elements, 
after detection of the i-th sector, are equal to the 
elements iL, . . . ,(i+l)L-l of i, as well as the vector 
d, which, after the detection of the i-th sector, is 



<h o {0h'0>- 

With the formal notation 30 



ai'.io=<'#;(i> 

the following is obtained 35 



0- 

The 4> hA (l p l'.i,i") depend only on the properties of the trans- 
mission channel. Trie way in which these can be determined 
in the receiver is described further below. 

The adaptation algorithm can now be summarized as 
follows. 

The received burst is split into 2N+1 sectors . . . 
each having M signal samples. It is assumed that the 
channel impulse response within a sector is constant. 

A vector £ of dimension (2N+1)L is defined, whose 
contents at any time are the estimated channel impulse 
responses of the individual sectors. The components 
£(iL), . . . ,x((i+l)L-l) of this vector are used as the 
channel impulse response for the detection of the i-th 
sector. Before starting the detection, a(iL+l) is initial- 
ized with h(l), the 1-th component of the channel 
impulse response estimated from the test sequence. 

A matrix E of dimension (2N+l)Lx(2N+l)L is defined, 
whose (iL+UlffO-th element is initialized before the 
start of detection using 



d=rrSfLi 

25 where i t is a vector of dimension M and contains the 
received signal samples associated with the i-th sector A,. 
After the detection of the i-th sector, 



is now calculated, and the recursions 

£ .^E+PSEt* are carried out 

These recursions result from trivial transformations of the 
general equations (8), (9) and (10), taking account of the 
model parameters. 

The particular advantage of this algorithm is justified by 
the fact that the channel impulse responses of all the sectors 
are updated in every recursion step. The k-th recursion of the 
channel impulse response in the sector A t in this case 
depends on the channel impulse responses of all the sectors 
in the recursion step k-1. This results in the capability of the 
algorithm to adapt to fading incursions. 

As stated above, in order to operate optimally, the algo- 
rithm needs to know the channel correlation properties 
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65 



However, it is found that the algorithm also operates satis- 
factorily if, instead of the above (^(U'iU'), completely 
different values, but values which are selected in a physi- 
cally sensible maimer, are selected for initialization of the 
error correlation matrix. $ M Q,r,i,i') depends in general on 
the Doppler frequency 



(that is to say, finally, on the vehicle speed and the carrier 
frequency) and the incidence angles of the received radiation 
components with respect to the movement direction of the 
mobile subscriber. Physically sensible initialization of the 
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error correlation matrix is obtained by estimating (^(U'fU') 
for the worst case situation and using these values for the 
initialization. 

Another option is to calculate, in the receiver, estimated 
values for the ^(l.l'.i.i"). For this purpose, the receiving 5 
process is started using the above worst case estimate for the 
$a/,(U\U'). After a burst has been completely detected, 
estimated values for the channel impulse responses h ; are 
available for all the sectors, and to be precise 2(iL+l), the 
(iL+l)-th component of the estimated state vector x (after 10 
the last step in the channel impulse response adaptation), is 
a suitable estimated value for h,(l). A suitable estimated 
value for <^ AA (l 1 l , ,i^ r ) is now obtained by averaging x(iL+ 
l)£*(i'L+r) over time. 
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I claim: 

1. A method for adaptation to the time-variable channel 40 
impulse response during a burst in receivers for TDMA 
mobile radio systems, each burst including in the center 
between two data regions a test sequence which is known to 
the receiver and by means of which an initial estimate of the 
channel impulse response is possible, wherein 45 
the received signal of a burst is sampled, 
the received burst is divided into 2N+1 sectors A_ N , ... 
,A +;v , in such a manner that the channel impulse 
response can be regarded as being constant within a 5Q 
sector, 

a vector X of dimension (2N+l)Lis defined, L being the 
number of channel impulse response coefficients and 
the (iL+l)-th component of X -N^i^+N, 0^1 iL-1 
is initialized before the start of detection using the 1-th 55 
coefficient ' of the estimated value of the channel 
impulse response obtained from the test sequence, 



a matrix E of dimension (2N+l)Lx(2N+l)L is defined, 
whose (iL+l,iI-+r)-th element is initialized using 



where ^ ll (l,r M ' I i,i"") is an approximation to <h t -(l)h/(l"")> 
and is the correlation of the 1-th component of the channel 
impulse response in the sector A, and the l'-th component of 
the channel impulse response in the sector A r , and a 2 /, is the 
variance of the error in the estimate of the channel impulse 
response coefficients from the test signal, 
a matrix E ; of dimension (2N+l)LxL is defined, whose 
columns, after the detection of the i-th sector, are the 
columns iL, . . . ,(i+l)LM of the matrix E, 
a matrix E a of dimension LxL is defined, whose lines, 
after the detection of the i-th sector are the lines iL, . . 
. ,(i+l)L-l of the matrix E„ 
a matrix S of dimension MxL is defined, whose (m,l)-th 
digit, after detection of the i-th sector, is occupied by 
reconstructed transmitted signal sample s(iM+m-l) 
where M is the number of signal samples in a section, 
a vector ^ of dimension L is defined, whose elements, 
after detection of the i-th sector, are equal to the 
elements iL, . . . .(i-t-l)L— 1 of x, 
a vector d of dimension M is defined which, after the 
detection of the i-th sector, is 



d=r r SH i 

where r,- is a vector of dimension M and contains the 
received signal samples associated with the i-th sector A;, 
each burst is processed in 2N steps, where 
in the k-th step, the i-th sector is detected using the 
components iL, . . . ,(i+l)L-l of the vector x, the sector 
A, being selected by a method by means of which, 
initially, all the data sections which are not affected by 
a fading incursion are provided for detection before an 
affected data section is processed, 
the matrix 



is calculated, where a 2 is the variance of the noise- contained 
in the received signal samples, 
and the recursions 



E ^e+pse; 
are carried out. 
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